To detect the looped Bloch bands of Bose-Einstein condensates in optical lattices 
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A loop structure was predicted to exist in the Bloch bands of Bose-Einstein condensates in optical 
lattices recently in [Wu and Niu, Phys. Rev. A 61, 023402 (2000)]. We discuss how to detect 
experimentally the looped band with an accelerating optical lattice through extensive and realistic 
numerical simulations. We find that the loop can be detected through observing either nonlinear 
Landau-Zener tunneling or destruction of Bloch oscillations. 
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I. INTRODUCTION 

The simple system of Bose-Einstein condensates 
(BECs) in optical lattices is of amazingly rich physics, 
as shown in recent theoretical and experimental studies. 
With certain choices of densities of BEC and strengths 
of optical lattice, this system exhibits various interesting 
phenomena, rangiiig from the dynamics of BEC Bloch 
waves ^ IE IE IS; Josephson effect 0, squeezed 
states Ig, and quantum phase transition between super- 
fluidity and Mott-insulator0. One can only expect more 
interesting physics to be discovered in this system, con- 
sidering the rich physics that we have known in the 
condensed-matter physics, where the prototype system 
is electrons in a crystal lattice. 

One very surprising finding for the system of BECs 
in optical lattices is the loop structure appeared in the 
Bloch bands as found in R.ef.[T(Tl|(see FigQ). This find- 
ing was later confirmed in further theoretical studies |E 
I 111 Il2| . This unusual and unique loop structure has 
very interesting physical consequences; First, it leads 
to the nonlinear Landau-Zener tunneling that is char- 
acteristically different from the linear Landau-Zener 
tunneling in particular, the nonzero tunneling in 
the adiabatic limit [Tof. Second, it destroys Bloch 
oscillationsHllllll. 

However, experimental exploration of this looped band 
and its related physical phenomena is yet to come forth. 
To our judgment, this experimental stalemate does not 
come from the lack of experimental techniques; it is the 
lack of enough theoretical guidance. On the one hand, 
it is not clear what signals to look for in an experiment 
to confirm the existence of the loop structure. On the 
other hand, one may be concerned that the unavoidable 
inhomogeneity of the BEC used in a real experiment may 
wash away all the interesting physics since the predicted 
loop structure and its related physics is based on the 
analysis of homogeneous BECs. 

The purpose of this paper is to dismiss these concerns. 
Based on extensive numerical calculations, we argue that 
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FIG. 1: Two lowest Bloch bands of a BEC in an optical lattice 
when C > V . The energy is in units of Ah^k\/m and the wave 
number k in units of 2fcL, where fcL is the wave number of the 
laser light that generates the optical lattice. 



one can confirm the existence of the loop structure in an 
experiment that involves dragging a high density BEC 
with an accelerating optical lattice. Similar experiments 
with low density BECs 5j or cold- atoms have been 
carried out to observe Bloch oscillations and Landau- 
Zener tunneling. As we will see later, there are two ways 
to look for the signs of the loop structure: the destruc- 
tion of Bloch oscillations and the observation of nonlin- 
ear Landau-Zener tunneling that shows a very distinct 
behavior from the well-known Landau-Zener tunneling. 

We focus on the experimental situations similar to 
what is described in Ref. [E S ^3 j where a realized BEC 
is loaded into an optical lattice and there is no trapping 
potential. These experiments can be well described by 
the mean-field Gross-Pitaevskii equation in one dimen- 
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where m is the atomic mass, fc^ is the wave vector of 
the laser hght that generates the optical lattice, Os is the 
s-wave scattering length between atoms, a is the accel- 
eration, and Vo is the strength of the potential which is 
proportional to the laser intensity. In our calculations, 
Gaussian functions are used to simulate the inhomoge- 
neous BECs loaded in optical lattices in the real exper- 
iments. Strictly, the lateral expansion of the BEG has 
certain effects on the longitudinal mot ion p^. In this pa- 
per we only consider the case where the lateral motion is 
negligible In experiments, another possible setup for 
quasi-one dimensional dynamics is to confine the lateral 
motion pi. fill . 

Without the acceleration, a = 0, the system becomes a 
nonlinear periodic system for which we can define Bloch 
waves as for a linear periodic system. 



(2) 



where ipk is periodic, ^pk{x + Tr/fci) — ^-)k{x). The Bloch 
state satisfies 



+VQC0s{2kLx)^k + 



IM'^^k- (3) 



The eigenenergies (or more precisely, the chemical poten- 
tials) fJ.{k) then form Bloch bands in the Brillouin zone. 
As shown in Fig^ when the interaction between atoms 
gets larger than a certain critical value, a loop structure 
is formed in the Bloch bands. 

In the following, we first briefly describe our numeri- 
cal methods. Then, we present our numerical results on 
nonlinear Landau-Zener tunneling and Bloch oscillations, 
and explain how to look for signs of the loop structure 
through these two phenomena. In the end, we discuss 
the relevance of our results to the experiments. 



II. NUMERICAL METHOD 

For the convenience of numerical calculations, we cast 
Eq.(^ into a dimensionless form 

where we have used the following set of scaled variables. 
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with no being the peak density of the BEG cloud. In 
writing down Eq.Q, we dropped the tildes (replacing 
X by X, etc.) without causing confusion. We use the 
Grank-Nicholson method for the numerical solution of 
Eq. (gl). This method preserves the unitarity of the time- 
evolution, and yields good convergence of the solutions 
for moderate values of the coupling strength C. Note that 
for the experiment where the lateral motion is confined, 
C has different definitions ^see one example in Ref . [TtI . 

In many experiments |j, l^, a BEG cloud has a typical 
size of order 10 /im, covering 100-200 wells of an opti- 
cal lattice. To model such inhomogeneous BEG clouds, 
we use a Gaussian wave packet as the initial state and 
then turn it into a inhomogeneous BEG Bloch wave by 
adiabatically turning on the optical lattice. 

The lattice strength V is taken to be smaller than 0.4 
(or 3.2 in units of the recoil energy 1i'^k\/2m). This 
choice serves two purposes. First, it guarantees that 
there be only one bound state inside each well so that 
we can use the separation of a BEG cloud to measure 
the tunneling probability as we will explain below. Sec- 
ond, it means that the mean field theory |p[J is a good 
description of the BEG system. 



III. NONLINEAR LANDAU-ZENER 
TUNNELING 

The nonlinear Landau-Zener tunneling has been stud- 
ied quite extensively in Ref.^^i where it is found 
to be ve_ry_ different from the linear Landau-Zener 
tunneling 13]. In particular, when C > V, the tunnel- 
ing probability no longer depends on the sweeping rate 
exponentially (sweeping rate is the acceleration for the 
system of a BEG in an accelerated lattice). Moreover, 
in the adiabatic limit where the acceleration approaches 
zero, the tunneling probability approaches a finite value, 
instead of zero as in the linear case. This non-zero tunnel- 
ing probability in the adiabatic limit is the direct result 
of the loop structure: when a Bloch state is driven to the 
edg e of the loop, it has to split, resulting in tunneling 
llOl . Therefore, experimental observation of this adia- 
batic tunneling, along with the non-exponential depen- 
dence of tunneling probability on the acceleration, can 
be viewed as a direct evidence of the looped band. 

However, the analysis in Ref. ^3 is based on a sim- 
plified two-level model derived with the assumption of 
homogeneity of BEG. This may leave experimentalists 
wonder whether the unique characteristics of nonlinear 
Landau-Zener tunneling can be observed in real experi- 
ments, where BEGs are inhomogeneous and span only a 
finite range of space. This concern is legitimate, but our 
numerical results show that the inhomogeneity does not 
blur up the essential physics. 

In our numerical simulations, we calculate the tunnel- 
ing probability through the separation of a BEG cloud. 
As we mentioned earlier, the lattice strength is chosen 
such that there is a only one band below the well barrier. 
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FIG. 2: Tunneling probability as a function of acceleration a 
for various values of C and V — 0.2. The case of C = is 
compared to the exponential function g-^-'^^^/'^ indicated by 
solid line without (•). 



As a result, the part of a wave packet tunneled into the 
upper band will not be dragged along the lattice while the 
part remained in the lower band will be dragged along. 
This leads to a separation of a BEC cloud after a certain 
time of acceleration. By integrating the left-behind wave 
packet, we obtain the tunneling probability. This tech- 
nique was actually used in experiments to measure the 
tunneling probability [5llT3|. 

Fig. |21shows our results of the nonlinear Landau-Zener 
tunneling probability with V — 0.2 for various values of 
C. The initial Gaussian wave packet is (j>{x,t = 0) = 

2/2 

, where uj is the width of the condensate and its 
typical values used is uj = 325. Optical lattice potential 
is turned on from V{x) = to V{x) — V for t < 40 
to achieve an inhomogeneous BEC Bloch wave. Then it 
is boosted with an acceleration a for two Bloch periods, 
and moves with a constant velocity afterward. 

In general, tunneling probability is greater for larger C 
at a fixed acceleration a, and greater for large accelera- 
tions a at a fixed C. For C < V, the tunneling probability 
changes with acceleration a rather exponentially and in 
particular, it goes to zero when a — *■ 0. This is exactly 
what is predicted in Ref. lOj. Therefore, this result im- 
plies that the inhomogeneity has only negligible impact 
on the tunneling probability when C < V. 

When C > V,we also do not see much impact of the in- 
homogeneity of BEC: the tunneling probability does not 
depend on the acceleration a exponentially and in partic- 
ular, when the data is extrapolated to the adiabatic limit, 
a — !■ 0, the tunneling probability seemingly tends to a 
non-zero value. This clearly shows that these two unique 



characteristics of nonlinear Landau-Zener tunneling are 
not washed away by the inhomogeneity of the BEC cloud. 
The experimental observation of them will be viewed as 
the confirmation of the loop structure. However, since 
the acceleration can not be made arbitrarily small in a 
real experiment and our simulations, one may argue the 
validity and confidence on the extrapolation of the data 
to the adiabatic limit of zero acceleration. We believe 
that this problem can be partially solved by repeating 
the experiment with many different densities of BECs. 
Consistency among different sets of data will be a strong 
support of the validity of extrapolation. On the other 
hand, the observation of non-exponential dependence of 
the tunneling probability on the acceleration will provide 
an unambiguous evidence. 



IV. BLOCH OSCILLATIONS 

Bloch oscillations occur when a Bloch state is driven 
across the Brillouin zone by a small external field (l^. 
For a BEC in an optical lattice, Bloch oscillations can 
be achieved by dragging the lattice with a small acceler- 
ation as reported experimentally in Refs.j^Q- In order 
to observe these oscillations, a key requirement is that the 
acceleration must be small enough so that the tunneling 
to the upper band is negligible. If tunneling probability 
into the upper band is increased, say, by large accelera- 
tions, Bloch oscillations can be destroyed. 

What is interesting with BECs is that we have another 
way to increase the tunneling probability besides increas- 
ing the acceleration. It is to increase the density of the 
BEC. Especially, as discussed in the last section, when 
the density is high enough such that C > V, there is a 
non-zero lower limit on the tunneling probability as the 
result of the loop structure in the energy band. In other 
words, no matter how small the acceleration is, we will 
not be able to observe Bloch oscillations due to the non- 
zero adiabatic tunneling when C is big enough. There- 
fore, the observation of breakdown of Bloch oscillations 
will provide another way to detect the loop structure. 

In an experiment, one can repeat the measurements 
with increasing densities of BECs for a fixed small accel- 
eration. One expects to observe Bloch oscillations when 
the density is low; as the density increases, the oscilla- 
tions will deteriorate and eventually be destroyed. This 
is indeed the case as shown in Fig. |2| from our simula- 
tions. In this figure, we show currents of the condensate, 
j = J {■^)lm{(f)* ^)dx, as a function of time in the accel- 
erating frame for a small a. For C < V, Bloch oscillations 
are preserved during the first Bloch period; however, for 
C > V, Bloch oscillations are disrupted and completely 
destroyed during the first Bloch period. Even for the case 
of C = Bloch oscillations are seriously disrupted. 

In our view, this is a much better way to detect the 
loop structure shown in Fig^than the method discussed 
in the last section regarding the nonlinear Landau-Zener 
tunneling. With increasing densities, the measurement 
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FIG. 3: Top panel: Currents as a function of time for various 
values of C at V = 0.0875 and a = 0.005. Bloch oscillations 
are destroyed clearly when C > V. Results are shown for 
C=0.02, 0.0875, and 0.15 in the accelerating frame. Bottom 
panel: The wave functions at different points of breakdown 
of Bloch oscillations, as indicated in the top panel at t = 
251.6(A), 283.0(B), 314.5(C), and 345.9(D). The smoothness 
of the wave functions shows no signs of dynamical instability, 
implying the destruction of Bloch oscillations is the result of 
nonlinear Landau-Zener tunneling. 



FIG. 4: Top panel: Currents as a function of time with V 
= 0.3175 and a = 0.02. Results arc shown for C = 0.2, 0.3, 
and 0.4 in the accelerating frame, and we see the destruc- 
tion of Bloch oscillations even when C < V, indicating that 
it is caused by the dynamical instability. Bottom panel: 
Densities of the wave function, |<^(a;)p, for C = 0.3 at 
t = 101.5(A), 118.4(B), 135.3(C), and 202.9(D). The "spiky" 
or "messy" wave functions at points C and D are the result 
of dynamical instability. 



the density the faster the cloud expands. This makes the 

of the tunneling probability based on the separation of a separation more difficult. With Bloch oscillations, we do 
BEC cloud will become more difficult because the higher not need to worry about this difficulty. 
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One caution we need to take is that there is an- 
other mechanism for breakdown of Bloch oscillations, 
the dynamical instability of Bloch waves discussed in 
Rcf. 2, 6, 12]. The instability will cause the system stray 
away from the Bloch states with very small amount of 
perturbations or noises. This is demonstrated in FigQJ 
where we see the destruction of Bloch oscillations even 
when C < V for a very small acceleration. The "spiky" 
or "messy" wavefunctions in the bottom panel signal the 
onset of dynamical instability. 

The dynamical instability always exists in the neigh- 
borhood of C ~ where the breakdown of Bloch 
oscillations starts to take place. As a result, one must 
make sure that the breakdown of Bloch oscillations is 
the result of the loop structure, instead of the dynamical 
instability. One way is to use weak optical lattices (small 
values of V). In this case, the growth rates of unstable 
modes are small therefore the dynamical instability will 
not dominate in the first few oscillations as seen in the 
bottom panel of Fig|31 These smooth wave functions 
indicate that the dynamical instability is yet to come 
into play. Therefore, the destruction of Bloch oscillations 
in the upper panel is purely due to the nonzero adiabatic 
tunneling resulted from the loop structure. 



the loop structure in the BEG Bloch bands. We sug- 
gested two possible scenarios: the observation of break- 
down of Bloch oscillations and nonlinear Landau-Zener 
tunneling in an accelerating lattice. Experiments sim- 
ilar to what we suggest have already been carried out 
with low density BECsQ, where the corresponding cou- 
pling strength is in the range of C = 0.026 ^ 0.04. 
The values used in our calculations cover the range of 
C — 0.05 ~ 0.3. This means that the signatures of loop 
structure studied in this paper can be observed by in- 
creasing the density of BECs by two to ten times. This is 
certainly possible with the current experimental set-ups. 
The density of the BEG in Ref.||], lO^cm-^, can be in- 
creased to meet the requirement considering the BEGs 
of densi ty a s high as 3 x lO'^^cm"'^ have already been 
achievedflSj. 
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